The statistical model in Pb-Pb collisions at the LHC 
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We briefly review the predictions of the thermal model for hadron production in comparison to 
latest data from RHIC and extrapolate the calculations to LHC energy. Our main emphasis is 
to confront the model predictions with the recently released data from ALICE at the LHC. This 
comparison reveals an apparent anomaly for protons and anti-protons which we discuss briefly. 
We also demonstrate that our statistical hadronization predictions for 3/iff production agree very 
well with the most recent LHC data, lending support to the picture in which there is complete 
charmonium melting in the quark-gluon plasma (QGP) followed by statistical generation of J/ifr 
mesons at the phase boundary. 



The quark-gluon plasma produced in ultra-relativistic nuclear collisions profoundly influ- 
ences the production of hadrons E[2]|4). For central collisions between large nuclei the yields 
of hadrons made up of light quarks can be described very well from AGS to RHIC energies 
(see and refs. therein) within the statistical model, with the chemical freeze-out temperature 
T, the baryo-chemical potential yU/, and the fireball volume V as the only parameters. Since the 
extracted temperature values, which first increase sharply with increasing beam energy, level 
off near r=160 MeV for energies ~\/s NN >20 GeV, while the baryochemical potential decreases 
smoothly as a function of energy, the extrapolation to LHC energy is straightforward. Conse- 
quently, analysis of the recently released hadron production data from ALICE at the LHC pro- 
vides a crucial test for the statistical model. This was first investigated in |5| for data on hadrons 
with light and heavy quarks. Here we provide an update and compare to the most recent data. 

Before proceeding to LHC data we show, in Fig. [TJ the situation for full RHIC energy 



( V^aw =200 GeV). The overall trend of the data is very well reproduced by the calculations. 
However, the yield of protons and anti-protons from PHENIX and Brahms is overpredicted, 
while the yield of multi-strange baryons is generally underpredicted, leading to a rather poor^ 2 
value of the fiQ We note that feeding from weak decays of multi-strange baryons could be much 
improved with vertex detectors, as have been (PHENIX) or will be (STAR) installed into the 
RHIC detectors. Measurements with vertex detectors would be very important to provide a more 
uniform experimental picture. 

In Fig. [2] we show the results of thermal model fits to reently released ALICE data [6|. The 
left panel shows a fit to all currently available data. The unexpectedly low yields for protons 



1 A fit to data from the STAR experiment alone leads to a significantly improved fit. 
Preprint submitted to Elsevier October 30, 2012 
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Figure 1: Comparison of thermal model predictions with RHIC data.The data are as compiled in |3|, with arecent update 
taking into account all available information on feeding via weak decays of multi-strange baryons. 



and anti-protons drive the temperature of the fit to a rather low value (T - 152 MeV) while the 
yield of multi-strange baryons is significantly underpredicted. This is somewhat similar to the 
situation observed at RHIC (Fig. [TJ. With the more than a factor of 2 smaller error bars of the 
ALICE data compared to results from the RHIC experiments the reduced^- 2 value approaches 4, 
and the temperature parameter is significantly lower than expected from the extrapolation from 
the data at lower energies 0. 

The right hand panel in Fig. [2] shows the result of excluding protons and anti -protons from 
the fit. This leads to a very good description of all remaining data, with excellent x 1 parameter 
and a temperature value (164 MeV) completely in line with expectations. Naturally, the nucleon 
yields are now about a factor of 1.4 below the calculated values. This apparent proton anomaly 
could be due to annihilation in the hadronic phase near the phase boundary. Indeed, schematic 
model calculations indicate such an effect We note, however, that annihilation affects 

not only nucleons, but also strange and multi-strange baryons. If annihilation is the explanation 
for the proton anomaly then the new ALICE data suggests that the annihilation rate for strange 
baryons is significantly less than that for nucleons. Further precision measurements, including 
also correlations among baryons and anti-baryons, are needed to shed light on this observation. 

In the following we use the statistical model to make predictions for charmonium production 
and compare the results to the most recent ALICE data ll9l lT0ll . Suppression of J /if/ mesons in 
the QGP was originally predicted ifTTIl as a key signature for a dense partonic phase. In contrast, 
in lfT2ll it was argued that charmonium production can be well described in the statistical model 
by assuming that all charm quarks are produced in initial, hard collisions. An important further 
input is that the QGP provides complete color screening, implying that charmed hadrons and 
charmonia are first produced at the phase boundary with statistical weights (for a recent review 
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Figure 2: Thermal model fits to ALICE data on hadron production in central Pb-Pb collisions. The left panel shows the 
result of the fit to all available data, while protons and anti-protons are excluded from the fit shown in the right panel. 
The ALICE data are preliminary results shown at this conference (3). 



see fOl . for a detailed more technical description see [14]). An important element is thermal 
equilibration of charm quarks, at least near the transition temperature T c , The new ALICE data 
|fT31 on spectra and flow of open charm hadrons and charmonia provide good evidence for this. 
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Figure 3: Centrality dependence of R J J^ for RHIC and LHC energies at mid-rapidity (left panel) and forward rapidity 
(right panel). The two curves shown for the LHC energy coorespond to a range of expected shadowing. The ALICE data 
shown in the left panel are preliminary results shown at this conference (9). 



The centrality dependence of the nuclear modification factor R J !f as measured recently by 
ALICE IflOl is shown in Fig. [3] for central and forward rapidity, and compared to RHIC data 
from the PHENIX collaboration |[T6l as well as to predictions from the statistical hadronization 
model. We first note that, at LHC energy, much less suppression is observed compared to the 
RHIC results, both at forward- and at mid-rapidity. The model calculations iTPfl reproduce this 



trend very well. In our model the larger Rjf values at midrapidity are due to the enhanced 
generation of charmonium around mid-rapidity, determined by the rapidity dependence of the 
charm production cross section. Also the observed centrality dependence is correctly reproduced. 

The successful description of the new ALICE data lends strong support to the interpretation 
that, at LHC energy, }/if/ mesons do not form or survive inside the QGP, implying strong color 
screening. Rather, the observations are consistent with the formation of charmonium bound 
states at hadronization of the QGP. Conceptually, this is very different from the mechanism of 
continuous formation and destruction of charmonia in the QGP, as employed in transport models 
ifTTl [181 . In our model, charmonium production is a direct signal for deconfinement of charm 
quarks: the charmonia are dominantly formed from initially uncorrelated c and c quarks. Further 
measurements of iff' and^ e states, as planned with the ALICE upgrade project lfl9ll will be crucial 
to differentiate between the models. The next step is a measurement in pPb collisions, which will 
clarify the contribution of shadowing. In Pb-Pb collisions measurements will be performed at 
full LHC energy. Due to the increase of the charm production cross section, we expect a further 
increase of R J JJ^ of up to 40% at mid-rapidity for central collisions l20l . 
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